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Abstract. We investigate the brightness evolution of 7 FU Orionis systems in the 1— 100 /xm wavelength range using data 
from the Infrared Space Observatory (ISO). The ISO measurements were supplemented with 2MASS and MSX observations 
performed in the same years as the ISO mission (1995—98). The spectral energy distributions (SEDs) based on these data 
points were compared with earlier ones derived from the IRAS photometry as well as from ground-based observations carried 
out around the epoch 1983. In 3 cases (ZCMa, Parsamian21, VI 331 Cyg) no difference between the two epochs was seen within 
the measurement uncertainties. V1057 Cyg, V1515 Cyg and V1735 Cyg have become fainter at near-infrared wavelengths while 
V346 Nor has become slightly brighter. V 1057 Cyg exhibits a similar flux change also in the mid-infrared. At A > 60 /im most 
of the sources remained constant; only V346Nor seems to fade. Our data on the long-term evolution of V1057Cyg agree 
with the model predictions of Kenyon & Hartmann (1991) and Turner et al. (1997) at near- and mid-infrared wavelengths, but 
disagree at A > 25 /im. We discuss if this observational result at far-infrared wavelengths could be understood in the framework 
of the existing models. 
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1. Introduction 

FU Orionis objects are low mass pre-main sequence stars de- 
fined as a class by Herbig (1977). Most members of this class 
have undergone outburst in optical light of 4 mag or more, fol- 
lowed by a fading phase on the timescale of several decades. 
In some cases (e.g. ZCMa, L1551 IRS 5, BBW 76) eruption in 
the optical was not observed but they were identified as class 
members on the basis of their spectral characteristics. For a 
comprehensive review of FU Orionis objects see Hartmann & 
Kenyon (1996). 

According to the most widely accepted picture the 
FU Orionis outburst is a consequence of a rapid temporal 
increase of the disk accretion rate (Hartmann & Kenyon 
1996). Predictions of these types of models (e.g. Kenyon & 
Hartmann 1991, Bell et al. 1995, Turner et al. 1997) have to 
be confronted with multiwavelength monitoring observations 
of the outburst period and the fading phase. The fading phase 
is well documented in the optical/near-infrared (Kopatskaya 
1984, Ibragimov & Shevchenko 1988, Kenyon & Hartmann 



1991). Very few data have been available so far at mid- and 
far-infrared wavelengths where thermal emission of the disk 
and the circumstellar envelope can be observed. Recently 
the Infrared Space Observatory (ISO, 1995—98, Kessler et 
al. 1996) provided new photometric data in the 4.8—200 /im 
range on FU Orionis-type stars. 
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In our study we search for systematic brightness variations 
during the post-outburst phases of 7 FU Orionis objects at mid- 
and far-IR wavelengths. We compiled two SEDs for each ob- 
ject: the first one is based on IRAS photometry as well as 
ground-based observations and is representative of the evolu- 
tionary status around 1983. The second SED was compiled 
from observations taken around 1996—2000 including data 
from ISOPHOT the infrared photometer on-board ISO (Lemke 
et al. 1996), MSX (Egan & Price 1996), and 2MASS (Cutri 
et al. 2003). The comparison of the two SEDs provides infor- 
mation on the wavelength dependence of the far-infrared flux 
evolution during a period of 15 years. As an example in the 
case of VI 057 Cyg we present a detailed comparison of our re- 
sults with the predictions of the models of Kenyon & Hartmann 
(1991) and Turner et al. (1997). 
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Object 


Outburst 


D[pc] 


IRAS 


MSX 


ISO 


2MASS 


BBW76 1 


<1930 ib 


1700 18 


07486-3258 


G248.7075-03.3686 




07503560-3306238 


CB34V 2 


<1994 2 


1500 18 




G186.9520-03.8325 




05470377+2100347 


FUOri 3 


1937 15 


450 18 


05426+0903 




SWS,LWS 


05452235+0904123 


V1057Cyg 3 


1970 15 


600 18 , 700 19 




G085.4595-01.0468 


PHT,LWS 


20585371+4415283 


V1515Cyg 3 


1950s 15 


1000 18 , 1050 3 


20220+4202 


G079.9187+02.7391 


PHT 


20234802+4212257 


L1551 IRS5 4 


? 


140 18 


04287+1801 






04313407+1808049 


PP13S 5 


<1900 16 


300 18 


04073+3800 






04104088+3807517 


Parsamian21 6 


? 


1800 18 ,400 20 


19266+0932 


G045.8149-03.8309 


PHT 


19290085+0938429 


Re50NIRSl 7 


1960-70 17 


460 18 






PHT,SWS 


? 


V1331Cyg 8 


7 


550 18 , 700 19 


20595+5009 


G090.3 121+02.6774 


LWS 


21010920+5021445 


V1735 Cyg 9 


1957-65 15 


900 18 


21454+4718 


G093.7587-04.6371 


PHT,LWS,SWS 


21472065+4732035 


V346Nor 10 


~1984 15 


700 18 


16289-4449 


G338.5458+02.1178 


CAM,PHT,LWS,SWS 


16323219-4455306 


V883 Ori 7 


? 


460 18 


05358-0704 






05381810-0702259 


ZCMa 11 


? 


930 18 , 1150 21 


07013-1128 


G224.6077-02.5574 


PHT,LWS,SWS 


07034316-1133062 


RNO IB 12 


? 


850 18 


00338+6312 


G121. 2940+00.6572 


CAM,LWS,SWS 


00364599+6328529 


RNO 1C 13 


? 


850 18 


00338+6312 


G121.2940+00.6572 


CAM,LWS,SWS 


00364659+6328574 


AR6A 22 


? 


800 22 




G203.2028+02.0653 




06405930+0935523 


AR6B 22 


? 


800 22 




G203.2028+02.0653 




06405930+0935523 


OO Ser 14 


1995 14 


31 1 23 






CAM,PHT,LWS,SWS 


18294913+0116206 


VI 647 Ori 24 


2004 


400 25 


05436-0007 




CAM 


05461313-0006048 



Table 1. Up-to-date catalogue of FU Orionis objects (for column description see Sect. l2.H . References: [1] Eisloffel et al. (1990); 
[2] candidate, other designation: VI 184 Tau, Yun et al. (1997); [3] Herbig (1977); [4] Can- et al. (1987); [5] Sandell & Aspin 
(1998); [6] Staude & Neckel (1992); [7] Reipurth 50 N IRS 1, Strom & Strom (1993), a 20 oo = 5 fe 40 m 27.4 s S 20 oo = -7°27'33"; 
[8] assumed to be in a pre-outburst state, McMuldroch et al. (1992); [9] Elias (1978); [10] Graham & Frogel (1985); [11] 
Hartmann et al. (1989); [12] Staude & Neckel (1991); [13] Kenyon et al. (1993); [14] Hodapp (1995); [15] Hartmann & Kenyon 
(1996); [16] Aspin & Sandell (2001); [17] Reipurth & Aspin (1997); [18] Sandell & Weintraub (2001); [19] Chavarria-K. (1981); 
[20] Eiroa & Hodapp (1990); [21] Herbst et al. (1978); [22] candidate, Aspin & Reipurth (2003); [23] de Lara et al. (1991); [24] 
FU Ori candidate, in the state of rapid brightening since January 2004. An associated submm source was detected by Mitchell et 
al. (2001); [25] Lis et al. (1999) 



2. Target list, observations and data reduction 

2. 1. Up-to-date catalogue of FU Orionis objects 

When Herbig (1977) defined the class of FUOrionis objects, 
there were only three known members. Since then a number of 
new objects have been identified. As a first step of our study we 
compiled an up-to-date list of FU Ori stars. We included FU Ori 
candidates without attempting to homogenise the identification 
criteria of different authors. The resulting catalogue, which is 
to a large extent similar to the list of Sandell & Weintraub 
(2001), is presented in Tab.^ In this table Col. 1 refers to the 
original paper where the object was identified as FU Ori type. 
The next two columns give the date of outburst (if known) and 
the distance to the object, respectively. Associated sources in 
the IRAS and MSX catalogues are given only in the case of 
good positional coincidence with the optical coordinates. The 
RNO 1B/1C system (separation ~5") was detected but unre- 
solved by both IRAS and MSX. Column 6 shows which of the 
four ISO instruments (ISOCAM, ISOPHOT, ISO-LWS, ISO- 
SWS) observed the source. These observations are available for 
the public in the ISO Data Archive 1 . The last column identifies 
the corresponding sources in the near-infrared 2MASS survey. 
These source names were derived from the sexagesimal coor- 



www.iso.vilspa.esa.es/ida/ 



dinates thus give the equatorial positions for the epoch 2000. 
(For the position of Re 50 N IRS 1 see figure caption.) 

2.2. Observations 

Table|2] lists the sources of infrared photometric data used in 
our study. The time distribution of the active periods of the 
instruments/projects provides a possibility to check for long- 
term variations of the infrared fluxes between ~ 1983 (IRAS, 
ground-based data) and 1996-2001 (ISO, MSX, 2MASS). 

From the sources in Tab.^we selected those 7 objects for 
further study where enough data were available at both epochs 
(-1983 and 1996-2000). Five of them have ISO multifilter 
observations to create a complete mid/far-infrared SED. One 
object, Parsamian21, was measured with ISOPHOT only at 65 
and 100 /im, but MSX data are available at shorter wavelengths. 
V 1 33 1 Cyg was not observed by ISOPHOT, but the availability 
of MSX fluxes makes possible the study of temporal variation 
at A < 25 [im. The detailed log of ISOPHOT observations is 
presented in Tab.0 

2.3. ISOPHOT data processing 

Most observations were carried out as pairs of ON and OFF 
measurements separated by 6—12'. The 180" aperture was 
selected because it was comparable to the IRAS far-infrared 
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Object 


Wavelengths [/im] 


Aperture ["] 


Date 


Obs. mode 


ISOJd 


V1057Cyg 


4.8,12,25,60,100 
65, 100 


180 
43x43 


17-Apr-1996 
8-Nov-1996 


ON/OFF 
PHT32 


15200607/15200608 
35800611 


V1515Cyg 


4.8, 12, 25, 60 
120 


180 
180x180 


24-Apr-1996 
24-Apr-1996 


ON/OFF 
ON/OFF 


15900605/15900606 
15900613/15900614 


Parsamian21 


65, 100 


43x43 


28-Sep-1996 


PHT32 


31601103 


V1735 Cyg 


4.8,12,25,60,100 
4.8, 12,25,60, 100 
65, 100 


180 
180 
43x43 


21-Apr-1996 
20-May-1996 
3-Dec-1996 


ON/OFF 
ON/OFF 
PHT32 


15600909/15600910 
18501409/18501410 
38301012 


V346 Nor 


4.8,12,25,60, 100 


180 


9-Feb-1996 


ON/OFF 


08402303/08402304 


ZCMa 


4.8, 12,25 
60, 100 


52 
99 


5-Nov-1997 
5-Nov-1997 


ON/OFF 
5x1 scan 


72003006/72003005 
72003007 



Table 3. Log of ISOPHOT observations. Observing modes are described in Sect. 12.31 ISO jd is the unique 8-digit identification 
number of ISO observations. 



Instrument 


Wavelengths [/im] 


Aperture 


Active period 


ground-based 


IHKLMNQ 


< 6" 


1970s- 


IRAS 


12, 25, 60, 100 


1-3' 


1983 


MSX 


4.25,4.29,8.28 


18" 


1996-1997 




12.13, 14.65,21.34 






2MASS 


I, H, Ks 




1997-2001 


ISOPHOT 


4.8-120 


43" -180" 


1995-1998 



Table 2. Sources of infrared photometric data used in our study. 



2003) we evaluated 23 observations of ISOPHOT standards 
(stars, asteroids, and planets) with default processing and de- 
rived an empirical correction curve from the comparison of the 
measured fluxes with their model predictions. The fluxes were 
extracted from the maps using an algorithm like the one de- 
veloped for ISOPHOT mini-map observing mode (Moor et al. 

2004) . The typical uncertainty of the empirical correction does 
not exceed the adopted 25% error bars. 



beam size. The far-infrared flux of Z CMa was extracted from 
a 5 x 1 scan obtained with the CI 00 camera. Three objects were 
measured at 65 and 100 in the high-resolution PHT32 map- 
ping mode (see The ISO Handbook vol. V, Laureijs et al. 2001). 

The data reduction was performed using the ISOPHOT 
Interactive Analysis Software Package VI 0.0 (PIA, Gabriel et 
al. 1997). After corrections for non-linearities of the integration 
ramps, the signals were transformed to a standard reset interval. 
Then an orbital dependent dark current was subtracted and cos- 
mic ray hits were removed. In case the signal did not fully sta- 
bilise during the measurement time due to detector transients, 
only the last part of the data stream was used. This was found 
mainly in observations with the 4.8 /im and 12 /im filters, while 
at other wavelengths the measurements showed sufficient sta- 
bility. The calibration of measurements at 4.8, 12, and 25 /im 
was performed by adopting the default responsivity. At longer 
wavelengths first we tried comparison with the on-board fine 
calibration source (FCS); however in most ON/OFF observa- 
tions the quality of the FCS measurement was unacceptable 
and we returned to the default responsivity. As error estimate 
we adopted an absolute calibration uncertainty of 25%, which 
represents well the sum of the random and systematic uncer- 
tainties. Colour corrections were applied for each measurement 
by convolving the observed SED with the ISOPHOT filter pro- 
file in an iterative way. 

Three objects (V1057Cyg, V1735Cyg, Parsamian21) 
were measured at 65 and 100 /im in the PHT32 observing mode 
(oversampled map). Since this method is not scientifically vali- 
dated we used our own procedure for evaluation. Instead of the 
highly interactive tool implemented in PIA (Tuffs & Gabriel 



2.4. IRAS data 

In order to have solid IRAS photometry for comparison with 
the ISOPHOT data we have re-analysed the IRAS raw data 
of two FUOrionis objects where the IRAS fluxes published 
by different authors (PSC, Weaver & Jones 1992, Kenyon & 
Hartmann 1991) were discrepant. We utilised the SCANPI pro- 
cessing tool at IPAC 2 . The results are listed in Tab.|4] remarks 
on individual objects are given below. 

V1 057 Cyg. This source was at the very edge of the IRAS 
coverage and only few scans are available. At 12—60 /im 
we adopted the fluxes published by Weaver & Jones (1992), 
though their formal uncertainties seem to be too optimistic. At 
100 /im, where only two detectors crossed the star, the source is 
surrounded by structures which are equally strong or stronger 
than VI 057 Cyg itself. The flux uncertainty is fully dominated 
by the accuracy of background determination. Experiencing 
with different background subtractions we estimate 47±10 Jy 
for V1057 Cyg's flux density at 100 /im. 

V1 51 5 Cyg. In this case the background is even more struc- 
tured than around VI 057 Cyg, and the coarse spatial resolution 
of IRAS is a problem already at 12 /im. At 12 and 25 /im we 
estimate 3.7±1.0Jy and 6.8±2.1 Jy, respectively, as integrated 
flux densities which also include an extended structure around 
the source, probably the infrared equivalent of the ring one can 
see in the optical wavelengths. The formal error bars reflect the 
uncertainty related to the background subtraction. At 60 /im 
the lack of spatial resolution is the limiting factor. The emis- 

2 http://irsa.ipac.caltech.edu/applications/Scanpi/ 
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Fig. 2. Spectral energy distribution of V1057Cyg. Open and 
filled symbols are as in Fig. ^ No reddening correction was 
performed. The lower panel shows the ratios of fluxes obtained 
at the two epochs (-1983 to 1996-2000). 



sion peak is offset by about 20" from the nominal position of 
V 1 5 1 5 Cyg, indicating that part of the emission from the region 
is unrelated to the star. For flux density towards the position of 
V1515Cyg our estimate is 25±10 Jy. At 100 /mi the emission 
is dominated by a local peak at about 1' from the position of 
V1515Cyg, and only an upper limit of 110 Jy can be derived 
for the FU Orionis object. 



3. Results 

The results of the ISOPHOT photometry as well as our newly 
derived IRAS flux densities are presented in Tab.0] Figure^ 
displays the SEDs of the sources from the two different epochs 
(-1983 and 1996-2000). The results for V1057Cyg, supple- 
mented with optical data, are shown in Fig.|2] 

By comparing the two SEDs of each source one can see the 
following evolutionary trends: 

- at near-IR wavelengths (A < 5 /im) the sources show var- 
ious behaviour. Parsamian21, VI 331 Cyg and ZCMa are 
unchanged, VI 057 Cyg, V 15 15 Cyg and VI 735 Cyg have 
faded, V346 Nor has become slightly brighter. 

- in the mid-IR (5 < A < 20 /im) only V1057 Cyg shows sys- 
tematic flux change: it became fainter by a factor of 2 dur- 
ing the period (lower panel in Fig.|2j. 

- at far-IR wavelengths (A > 60 /mi) five stars (V1057Cyg, 
VI 735 Cyg, ZCMa, Par 21, and V 15 15 Cyg) remained 
constant while V346 Nor seems to have become fainter. For 
VI 331 Cyg there are no FIR data other than IRAS. 

Below we make remarks on the individual sources. The case of 
VI 057 Cyg will be discussed in details in Sect.01 



V1 51 5 Cyg. According to the 850 /im submillimetre map of 
Sandell & Weintraub (2001) this source is situated in the mid- 
dle of a cold arcminute-size core, and this extended emis- 
sion is also visible on the MSX 8.28 /im map. In addition 
the 2MASS and MSX data reveal a compact object 2MASS 
20235198+4211260 at about 80" from V1515Cyg. Since this 
source is not included in the MSX Point Source Catalog, we 
extracted fluxes from the MSX maps using aperture photome- 
try. The resulting SED peaks at «4 /im and its position on the 
J— H vs. H— K diagram suggests a reddened background star 
with Ay > 15 mag. The contribution of this source in the 3' 
ISOPHOT beam can explain why the ISOPHOT 4.8 /im point is 
higher than the M-band observation of 1983. In the 12—25 /mi 
range ISO fluxes are higher while the MSX fluxes are lower 
than IRAS. This result, however, can be explained by the pres- 
ence of the arcminute core associated to V 15 15 Cyg, taking 
into account the different beam sizes of the three instruments 
(MSX: 18", IRAS: 1', ISOPHOT: 3'). At near-infrared wave- 
lengths the data indicate a slight (~ 15%) flux decrease, which 
is also visible in the K-band lightcurve (Fig.[3}. In the mid- 
infrared the comparison of the 1983 and 1996—98 SEDs are 
complicated by the beam confusion. The fact, however, that at 
—20 /im Q-band measurements from 1974 (4.1 ±0.7 Jy, Cohen 
1974) and 1989 (2.9±0.6 Jy, Kenyon et al. 1991) agree with the 
MSX flux value of 3.0±0.7Jy, suggests that the mid-infrared 
fluxes of V 15 15 Cyg remained constant. At far-infrared wave- 
lengths the emission of V 15 15 Cyg was constant within the 
measurement uncertainties. 

V1 735 Cyg. The 850 /im submillimetre map of Sandell & 
Weintraub (2001) reveals a cold compact source, V1735Cyg 
SM 1 at 20" from the star dominating the emission at both 450 
and 850 /im. From a careful analysis of the MSX images we 
concluded that in the whole near- and mid-infrared regime the 
submillimetre source is invisible and the observed emission can 
be assigned to V1735 Cyg. Thus, the star is responsible for the 
observed flux decrease of » 40% below 5 /zm, in accordance 
with the change in the K-band lightcurve (Fig.[3J. No flux vari- 
ation is visible at mid-infrared wavelengths. The relative con- 
tribution of the submm source, however, is probably increasing 
at longer wavelengths, as indicated by the fact that the IRAS 
position is located between the positions of the FU Orionis ob- 
ject and that of SM 1. Both the IRAS and ISOPHOT beams 
include both sources. Our result of no flux variation at far- 
infrared wavelengths indicates that the emission of V1735 Cyg 
was constant. 

V346 Nor. This star is situated in a complex region including 
the PMS star Reipurth 13, HH 56/57 and some cold extended 
emission within a radius of 60" (Sandell & Weintraub, 2001). 
Applying a maximum entropy based deconvolution technique 
Prusti et al. (1993) were able to determine the relative contribu- 
tion of Re 13 (separated by rj 50" from the FU Orionis star) in 
the IRAS bands, and found that its flux is a factor of 3-5 times 
lower than the emission of V346 Nor. At near-infrared wave- 
lengths the FU Orionis star - in a unique way in our sample - 
became brighter. This behaviour, reported already by Prusti et 
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Fig. 1. Spectral energy distributions of 6 FU Orionis-type objects. Open symbols correspond to an epoch of ~1983 while filled 
symbols mark the SED as of 1996—2000. The data are presented with no reddening correction. Error bars smaller than the symbol 
size were not plotted. For V1331 Cyg there are no NIR data from 1983; we plotted data from 1979 (the K-band lightcurve shows 
constant flux in this period). In the case of V346Nor the ISO and IRAS beams also include the nearby pre-main sequence star 
Re 13. For comparison the plotted MSX flux densities are sums of the individual fluxes of the two objects. 



al. (1993), is also shown by the K-band lightcurve (Fig.[3}. The 
brightening seems to be wavelength dependent: its amplitude 
increases monotonically from J to K. The same trend can still 
be observed at mid-infrared wavelengths up to 25 where the 
IRAS and ISOPHOT fluxes are the same. On the other hand at 
A > 60 /im the observed fluxes became lower between 1983 and 
1996. Assuming that Re 13 (which is also included in the beam) 
is not variable we may conclude that V346 Nor faded at 60 and 
100 fim. It is remarkable that the flux change is associated with 
a change in the far-infrared colour, too, which became warmer 
during the 15 years. The observed change could be an indica- 
tion for a temperature rise in the outer parts of the circumstellar 



environment, but it might also be an artifact caused by the dif- 
ferent sensitivities of the IRAS and ISOPHOT instruments for 
the extended emission. 



Z CMa is a binary system with a separation of 0. 1". The com- 
ponent which dominates at infrared wavelengths is probably 
a Herbig Ae/Be star, while its optically brighter companion is 
an FUOrionis object. The system is embedded in an elongated 
disk-like submm core (Sandell & Weintraub 2001), and was 
unresolved for both IRAS and ISOPHOT. According to the K- 
band lightcurve (Fig.|5J ZCMa did not change after 1974. The 
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Fig. 3. K-band lightcurve of 6 FUOrionis objects for the 
period 1963—2003. Data are taken from the literature (for 
Parsamian21 no observations are available). Grey stripes mark 
the two epochs (1983 and 1996—2000) which form the baseline 
for the long-term evolution. 



same conclusion can be drawn for the whole 2—100 /im SED. 
This result, however, does not exclude small flux variation of 
the FU Orionis component since the infrared SED of the system 
is dominated by the intermediate mass component. 

Parsamian 21 . This source is relatively isolated, no extended 
emission is associated with it at submillimetre wavelengths 
(Henning et al. 1998), thus beam differences of the different 
instruments are unimportant. The comparison of the two SEDs 
does not reveal any long term flux variation between 1983 and 
1996—98 within the measurement uncertainties. 

V1 331 Cyg is believed to be in a state prior to an FU Orionis- 
type outburst (Welin 1976, McMuldroch et al. 1993). The 
source is compact at submillimetre wavelengths (Henning et 
al. 1998), thus the comparison of the various observations is 
not influenced by beam effects. Our data show no evidence of 
temporal variation of the source's flux in the whole infrared 



regime. The K-band lightcurve (Fig.fJJ is consistent with this 
conclusion. 

4. Discussion: the case of V1 057 Cyg 

In this section we compare our new observational results on 
the temporal evolution of FU Orionis objects with predictions 
of models developed to describe their circumstellar structure. 
From our sample V1057Cyg has the best documented multi- 
wavelength flux evolution following its outburst in 1970. The 
ISO data points, supplemented with other infrared observations 
(Fig.fJJ, reveal a fast temporal evolution between ^1983 and 
1996—2000 (the fastest rate in the sample). The wavelength 
dependence of the evolution is also well determined. Since de- 
tailed models of VI 057 Cyg are available in the literature, we 
chose this object to compare the new observational results on 
the temporal evolution with model predictions. 

As far as we know there are two models of the circumstel- 
lar environment of VI 057 Cyg to fit the complete 1 — 100/im 
infrared SED (Kenyon & Hartmann 1991, Turner et al. 1997). 
Kenyon & Hartmann (1991) collected and analysed all data 
available until 1991 and modelled the post-outburst evolution 
in the framework of an accretion disk model. They concluded 
that a simple accretion disk model cannot explain the large far- 
IR fluxes, and investigated two possible geometries, a flared 
disk and an extended infalling envelope. Since the flared disk 
would require very high coverage factor, the authors favour a 
model where a flat disk is embedded in a spherically symmetric 
envelope with a wind-driven polar hole. The envelope, which is 
the remnant of the molecular cloud core, reprocesses the radi- 
ation from regions close to the central star. Turner et al. (1997) 
fitted the SED of VI 057 Cyg by computing outbursting flared 
disk models in which the mass flux varies with radius. The 
model includes reprocessing of disk emission by other parts 
of the disk. An envelope with a central hole exposing the inner 
disk is included as a layer of uniform thickness on the top of 
the disk. 

At optical and near-IR wavelengths (\<2.2 fim) both 
models claim that emission of the central source (the star plus 
the innermost part of the accretion disk) dominates the ob- 
served flux. After the outburst the accretion rate close to the star 
decreases leading to a flux decrease in this wavelength range. 
Our data confirm that the flux decay was observable between 
1983 and 2000, too. The data reveal that in the R, J, H, and K 
bands the flux dropped by the same constant factor of 2, while 
at B and V a larger decay was observed (lower panel of Fig.fJJ. 
This wavelength dependence may reflect the drop of effective 
temperature leading to a shift in the peak of the emission of the 
central source towards longer wavelengths. 

Between 3 and 10 /im the origin of the emission is the re- 
lease of accretion energy in the disk (Kenyon & Hartmann 
1991), and also starlight reprocessed in the same part of the 
disk and in an envelope (Turner et al. 1997). Decreasing accre- 
tion rate at the center (which leads to a decrease in the bolo- 
metric luminosity) would cause the drop of emission of all 
three components. The timescale at which viscous accretion 
could change is related to the dynamical timescale at this dis- 
tance (several years at ~1 AU, Pringle 1981). The timescale 
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A [/im] 


V1057Cyg 


V1515Cyg 


V1735Cyg 


V346Nor 


ZCMa 


Parsamian21 


V1331Cyg 


4.8 


2.67 ±0.67 


1.57 ±0.39 


1.29 ±0.32 


4.59 ± 1.15 


53.16 ± 13.29 






12 


5.68 ± 1.42 


11.64 ±2.91 


1.62 ±0.40 


6.61 ± 1.65 


97.33 ± 24.33 






25 


23.19 ±5.80 


11.01 ±2.75 


4.94 ± 1.24 


30.82 ± 7.71 


182.8 ±45.7 






60 


52.99 ± 14.11 


25.80 ±6.45 


41.76 ± 10.44 


46.54 ± 13.45 


290.1 ±72.5 






65* 


42.26 ± 10.57 




38.68 ± 9.67 






18.51 ±4.63 




100 


34.50 ±9.28 




77.23 ± 19.31 


39.52 ± 11.65 


479.1 ± 119.8 






100* 


62.06 ± 15.52 




86.82 ± 21.71 






20.71 ±5.18 




120 




78.45 ± 24.03 












IRAS- 12 


14.89 ±0.11 


3.7 ±1.0 


2.19 ±0.20 


9.73 ± 0.39 


125.1 ±5.0 


0.80 ±0.06 


1.12 ±0.03 


IRAS-25 


28.72 ±0.08 


6.8 ±2.1 


8.09 ±0.52 


30.98 ± 1.33 


204.5 ± 8.7 


4.07 ±0.26 


2.62 ±0.02 


IRAS-60 


53.71 ± 2.48 


25 ± 10 


40.84 ± 4.90 


69.07 ± 4.83 


312.3 ±40.6 


11.49 ± 1.03 


6.88 ±0.23 


IRAS- 100 


47.0 ± 10.0 


< 110 


92.97 ± 15.63 


74.91 ± 4.20 


375.5 ±63.2 


15.58 ± 2.04 


8.22 ± 1.19 



Table 4. Photometric results of ISOPHOT and IRAS. All fluxes are colour-corrected and presented in Jy units. (*): PHT32 ob- 
serving mode. Source of IRAS data: PSC - V1735Cyg, V346 Nor, Z CMa, Parsamian21; Weaver & Jones (1992) - VI 331 Cyg; 
this work (Sect.l2.4> - V1515Cyg. In the case of V1057Cyg the 12-60 /xm fluxes are from Weaver & Jones (1992), the flux at 
100 /im was re-determined in this work (Sect. l2.4t . 



on which the reprocessed emission follows the fading of the 
central source depends on the optical depth of the reprocess- 
ing medium. For optically thin medium (disk surface layer, en- 
velope) an almost instantaneous reaction is expected (several 
hours, Chiang & Goldreich 1997, Eq. 25). For optically thick 
emission (disk interior) the thermal timescale related to the gas 
component dominates (several years at ^1 AU for a typical T 
Tauri star, Chiang & Goldreich 1997, Eq. 28). In the case of 
V1057Cyg the 3 — 10/im emission is probably not completely 
optically thick, as indicated by the presence of weak spec- 
tral features (silicate emission at 9.7 /im, Wooden et al. 1995, 
Hanner et al. 1998). Since our temporal baseline (1983-1998) 
significantly exceeds the mentioned timescales and since the 
post-outburst evolution started already in the seventies, the 
model prediction is that the 3 — 10/im emission is decreasing 
synchronised with the rate of the optical/near-IR decay in a 
wavelength independent way. Our measurements are fully con- 
sistent with these predictions. (Fig.[2]lower panel) 

The emission at A > 10 /im is reprocessed starlight in both 
the model of Kenyon & Hartmann (1991) and of Turner et 
al. (1997). The infrared radiation emerges from an envelope 
where dust particles are in radiative equilibrium with the il- 
lumination from the central region (in the model of Turner 
et al. the envelope is optically thick for its own radiation at 
A < 30 /im). The temporal evolution implied by the models is 
similar to the predicted trend at shorter wavelengths: a wave- 
length independent fading at the same rate as the fading of the 
central illuminating source (although the time-lags are gradu- 
ally longer at longer wavelengths, Chiang & Goldreich 1997). 
Our results are not fully consistent with this picture. Though 
at 10/im the flux dropped between 1983 and 1996 by a fac- 
tor of 2 similarly to the optical/near-IR rate, at far-IR wave- 
lengths (A > 60 /im) the flux of V1057 Cyg remained constant. 
The comparison of the IRAS 60 and 100 /im photometric points 
with our new ISOPHOT data at the same wavelengths clearly 
demonstrates that the far-infrared fluxes of VI 057 Cyg showed 
no variation between 1983 and 1997. The transition between 
the two types of behaviour is around 25 /im. 



Our results indicate that - unlike in the above mentioned 
models - at A > 10 /im two important emission components 
have to be taken into consideration. Between 10 and 25 /im we 
probably observe the envelope, while at longer wavelengths the 
emission cannot be explained this way. The material responsi- 
ble for the A > 25 /im emission is apparently not an optically 
thin medium reprocessing the light of the central source, be- 
cause its radiation does not react immediately on the source's 
fading. Moreover, the emitting material has the following prop- 
erties: (1) it is relatively cold; (2) it has a flat spectrum below 
100 /im, suggesting aT~ r -°- 5 radial temperature profile. 

The first possibility for the nature of the far-infrared emit- 
ting component could be a reprocessing medium which is opti- 
cally thick (e.g. a flared disk), and its time-lag behind the fading 
of the central source significantly exceeds our temporal base- 
line. Kenyon & Hartmann (1991), however, demonstrated that 
the outer part of a flared disk around V 1 057 Cyg would produce 
an unacceptably large covering factor and rejected this possi- 
bility. In addition, based on more detailed calculations, Bell et 
al. (1997) pointed out that such a flared structure would not be 
realised in nature. At a radius of ~ 10 AU the flared disk would 
develop into a 'decreasing curvature' disk, corresponding to the 
condensation of water ice, thus the outer disk is shadowed from 
the illumination of the central source. Such a disk cannot pro- 
duce either the T ~ r ' 5 temperature profile or the high level 
of far-infrared emission. Based on these arguments we exclude 
reprocessing of the radiation of the central source in the outer 
part of the disk as the origin of the far-infrared emission. 

Nevertheless, it may be interesting to mention that the crit- 
ical radius where the disk turns to concave (i.e. the border line 
between the illuminated and shadowed parts of the disk) cor- 
responds to a temperature of «300K (Bell et al. 1997). Since 
the blackbody emission of this temperature contributes mainly 
to the 25 /im range of the SED, one might speculate that the 
change in the temporal behaviour of the SED of V1057Cyg, 
also at 25 /im, could be related to the turnover of the disk pro- 
file. However, in the previous paragraphs we concluded that 
neither the short wavelength emission (A<25/itn, attributed 
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an envelope) nor the long wavelength part of the SED are 
dominated by the radiation of an optically thick passive disk. 
We think that the agreement of the transition wavelength in 
VI 057 Cyg at 25 fim with the peak of the blackbody curve at 
the sublimation temperature of water ice in the disk is a ran- 
dom coincidence. This conclusion is further supported by the 
case of V1735 Cyg, where the shortest wavelength of constant 
flux is at 10 ^m, clearly different from the temperature of the 
turnover radius (Fig.Q. 

As a second possibility for the origin of the far-infrared 
emission Bell et al. (1997) proposed that the heating of the 
outer part of an optically thick disk could be dominated by the 
ambient radiation field of the cloud core around the young star, 
rather than by the radiation of the central source. As a result, 
the radial temperature profile becomes constant at the ambient 
temperature value, producing a shallower - or even approxi- 
mately flat - far-infrared SED. Since its origin is unrelated to 
the central source, the far-infrared emission is not expected to 
vary in time. 

A third possibility would be if the far-IR emission is pow- 
ered by accretion in the outer disk (active disk), where the 
accretion rate is constant and unrelated to the outburst of the 
central region (Bell et al. 1995). However, the spectral shape 
of such a standard accretion disk would follow the canonical 
vF v ~ A -4 / 3 law, in contradiction with the observed flat spec- 
trum (Fig.|2j. A possible solution for this problem could be if 
the temperature profile of the accretion disk differs from the 
standard one. Lodato & Bertin (2001) suggested that in a self- 
gravitating protostellar disk the non-Keplerian rotation curve 
may result in a nonstandard temperature distribution, and thus 
in a flat spectrum in the far-IR. The high disk mass required 
by this model (« 1 Mq) is not inconsistent with measurements 
of some FUOrionis objects (Henning et al. 1998, Sandell & 
Weintraub 2001), but the measured disk mass of VI 057 Cyg 
(« 0.1 Mq, Sandell & Weintraub 2001) seems to be too low 
causing any observable effect. 

The last possibility would be to assume that the origin of 
the far-IR emission is unrelated to VI 057 Cyg. Such a source 
could be a - so far undetected - embedded IR companion. It 
is interesting to note that Bell et al. (1995) already speculated 
about such a companion which could have triggered the out- 
burst of VI 057 Cyg by a close passage. A deeply embedded in- 
frared companion, like e.g. the very red object discovered in the 
vicinity of LkHa 198 at 10/^m (Lagage et al. 1993), might be 
too faint to be detected at near-infrared wavelengths so far. At 
mid- and far-infrared wavelengths VI 057 Cyg has never been 
imaged with high-resolution instruments, thus the existence of 
a companion unresolved in the beam of the MSX (18") cannot 
be excluded. 

In the previous paragraphs we argued that an optically thick 
disk, reprocessing the radiation of the central source, or accre- 
tion with a non-Keplerian rotation curve are less likely expla- 
nations, but a shallower temperature profile dictated by the am- 
bient radiation field or the possibility of an unrelated infrared 
companion cannot be excluded. The observations (Fig. 1) sug- 
gest that the constancy of far-infrared emission might be a gen- 
eral property of FU Orionis objects. Though a growing num- 
ber of FU Orionis objects are found to be binaries (Reipurth 



& Aspin 2004), it is not proved yet that all FU Orionis objects 
have companions. Thus we think that the most likely origin of 
the far-infrared radiation in VI 057 Cyg (and perhaps in other 
FU Orionis objects) is emission from the outer part of an opti- 
cally thick disk whose temperature profile is controlled by the 
ambient radiation field. 

Our investigation demonstrates that temporal variability of 
the infrared emission can provide important information on the 
physical properties of circumstellar material around FU Orionis 
objects. This kind of diagnostics of the disk structure - pro- 
posed first by Chiang & Goldreich (1997) - can also be applied 
on any other types of young stellar objects where the energy 
source of the disk, either accretion or reprocessing, is time- 
variable (e.g. exors). New infrared data on some FUOrionis 
objects are expected from the Spitzer Space Telescope which 
could show whether the evolutionary trends outlined by the ISO 
data also continue in our days, and whether the variation of the 
SEDs at longer wavelengths are responding to the flux drop at 
shorter wavelengths. 
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